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ABSTRACT 


We present radio observations of the afterglow of the bright 7 -ray burst 
GRB 991208 at frequencies of 1.4, 4.9 and 8.5 GHz, taken between two weeks 
and 300 days after the burst. The well-sampled radio light curve at 8.5 GHz 
shows that the peak flux density peaked about 10 days after the burst and 
decayed thereafter as a power-law Fr oc t-i- 07 ±o.o 9 ^ This decay rate is more 
shallow than the optical afterglow of GRB 991208 with Fq oc which 
was measured during the first week. These late-time data are combined with 
extensive optical, millimeter and centimeter measurements and htted to the 
standard relativistic blast wave model. In agreement with previous hndings, we 
hnd that an isotropic explosion in a constant density or wind-blown medium 
cannot explain these broadband data without modifying the assumption of a 
single power-law slope for the electron energy distribution. A jet-like expansion 
provides a reasonable £t to the data. In this case, the flatter radio light curve 
compared to the optical may be due to emission from an underlying host galaxy, 
or due to the blastwave making a transition to non-relativistic expansion. The 
model that best represents the data is a free-form model in which it is assumed 
that the broadband emission originates from a synchrotron spectrum, while the 
time-evolution of the break frequencies and peak flux density are solved for 
explicitly. Although the decay indices for most of the synchrotron parameters 
are similar to the jet model, the evolution of the cooling break is unusually rapid 
Uc oc t~‘^, and therefore requires some non-standard evolution in the shock. 

Subject headings: gamma rays: bursts - radio continuum: general - cosmology: 
observations 
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1. Introduction 

The afterglow of GRB 991208 was one of the brightest observed to date. Shortly after 
the coordinates of the 7 -ray error box were made known, a previously uncataloged radio 
source was discovered near its center. Based on the source location, its inverted spectrum 
and compactness, it was a strong candidate for being the radio counterpart of GRB 991208 
(Frail et al. 1999). This was conhrmed by the subsequent discovery of an optical transient 
at the same location as the radio source (Gastro-Tirado et al. 1999). Follow-up optical 
observations taken from 2 to 7 days after the burst showed that the spectral and temporal 
characteristics of this burst could be represented by power-law indices a = 2.2 ± 0 . 1 , and 
f3 = 0.75 ± 0.03, where F{t) oc (Sagar et al. 2000). With added data, Gastro-Tirado 

et al. (2001) determined a = 2.3 ± 0.07, and a somewhat steeper (3 = 1.05 ± 0.09. The 
7 -ray properties, the localization and the subsequent afterglow discovery for GRB 991208 
are presented in more detail by Hurley et al. (2000). 

GRB 991208 was also very bright at millimeter wavelengths (Shepherd et al. 1999) 
and as a consequence, well-sampled spectra and light curves covering frequencies between 
1.43 to 240 GHz were taken during the initial two weeks following the burst (Galama et 
al. 2000; Paper I). These observations made it possible, for the hrst time, to trace the 
time-evolution of the synchrotron spectrum directly without resorting to specihc models 
for the evolution of the characteristic frequencies and the peak flux density. In comparing 
these results with model predictions it was found that spherically symmetric explosions 
in homogeneous or wind-blown circumburst media could be ruled out. A model in which 
the relativistic outflow is collimated {i.e., a jet) was shown to account for the observed 
evolution of the synchrotron parameters, the rapid decay at optical wavelengths, and the 
observed radio to optical spectral flux distributions, provided that the jet transition has not 
been fully completed in the hrst two weeks after the event. 
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This high quality dataset stimulated other modeling efforts for GRB 991208. Li and 
Chevalier (2001) argued that an isotropic explosion in a wind-blown ambient medium 
described the data equally well. However, in order to account for the fast decay of the 
optical light curves with respect to the radio light curves, they invoked a non-standard 
break in the electron energy distribution. Panaitescu & Kumar (2002) also required a 
broken power-law and they concluded that the best £t to the data was obtained for a 
jet-like outflow with an opening angle of about 15° expanding into either a constant density 
or wind-blown circumburst medium. Dai & Gou (2001) present a variation on these other 
efforts, arguing that the steep optical decay is due to an anisotropic jet viewed off axis, 
which is expanding into a wind from a red supergiant. 

In view of the disparate conclusions reached on the basis of this early afterglow data, 
there is some hope that with further monitoring of the afterglow it may be possible to 
distinguish between these models. Accordingly, in this paper we present radio observations 
of GRB 991208 made with the VLA and the Very Large Baseline Array (VLBA) which 
began at the end of the two week period covered in Paper I and continued until the burst 
faded below our detection limits. We focus on the interpretation of the observations in 
terms of the predictions made by relativistic blast wave models. The observations and 
results are presented in §2 and §3 and are discussed in §4. 
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2. Observations 
2.1. Very Large Array 

Observations at 1.43, 4.86 and 8.46 GHz were made using the NRAO Very Large Array 
(VLA)^. All observations were performed in the default continuum mode in which, at each 
frequency, the full 100 MHz bandwidth was used in two adjacent 50-MHz bands. The 
synthesized beam at 8.46 GHz was approximately 0.8". The beam size at other frequencies 
scales inversely with frequency. Gain calibration was carried out by observing 3G48. The 
array phase was monitored by switching between the GRB and the phase calibrators 
J1637+462 (at 8.46 GHz), J1658+476 (at 4.86 GHz) and J1653+397 (at 1.43 GHz). Data 
calibration and imaging were carried out with the AIPS software package following standard 
practice. We performed phase self-calibration on the 1.43 GHz data. Flux densities were 
determined by a least squares htting of the Gaussian-shaped synthesized beam to the 
source. For very faint detection levels (< 3(j) we used the peak flux density at the position 
of the radio counterpart (Frail et al. 1999). A log of the observations is provided in Table 1. 


2.2. Very Long Baseline Array 

Observations centered at 8.35 and 4.85 GHz were made using the NRAO Very Long 
Baseline Array (VLBA) at four epochs. At both frequencies the observations consisted of 
eight channels of bandwidth 8 MHz (64 MHz total bandwidth) spaced equally over a ~ 0.5 
GHz frequency span. The phases were calibrated using the calibrator source 1637-1-472. 

At each epoch the observations lasted ~ 9 hours (divided over the source, calibrators 

^The NRAO is a facility of the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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and frequency switches). These observations were designed for an interstellar scintillation 
experiment (requiring substantial time and frequency span), the results of which will be 
presented in a separate paper. Here, in Table 2, we report just the time-averaged flux 
densities. 


3. Results 

Multiwavelength light curves of the event are presented in Figure 1. Data are from 
this work (Tables 1 and 2) and Paper I. Optical observations are taken from Castro-Tirado 
et al. (2001), were we have only included early data (less than 10 days after the event) to 
measure the ‘pure’ afterglow emission [i.e., where the underlying host does not contribute 
signihcantly). 

The VLBA measurements may suffer from amplitude losses as a result of imperfect 
phase referencing. The lightcurves indicate that such losses may be signihcant since at 8.46 
GHz the VLBA data seem to be systematically lower than the VLA data. 

The 8.46 and 4.86 GHz light curves rise to a maximum around 10 days. While the 
8.46 GHz light curve is smoothly rising to a maximum of about 2 mJy during the hrst 
week, it undergoes a sudden drop to about 1 mJy (around 9 days after the burst) followed 
by a second peak (around 12 days). Such erratic flux density variations are the hallmark 
of interstellar scintillation (ISS). We will return to this point in §3.1. The overall trend in 
Figure 1 is for the peak flux density to decline with decreasing frequency. A linear least 
squares £t was carried out to the late time radio light curve (At=53-293 days) at 8.46 GHz, 
and we derived a temporal decay index or = 1.07 ± 0.09 (y^^ = 1.07, where Fr oc 
This is to be compared with the much steeper optical decay cto = 2.2 ±0.1 determined at 
At=2-10 days (Sagar et al. 2000). Reconciling this difference between the optical and radio 
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decay indices is one of the main challenges in modeling the afterglow of GRB 991208 (see 

§4). 


3.1. Interstellar Scintillation 

The lines-of-sight towards most GRBs traverse a considerable pathlength through the 
turbulent ionized gas of our Galaxy. As a result, the radio emission from GRB afterglows is 
expected to be affected by interstellar scattering (ISS). Indeed, as predicted by Goodman 
(1997), ISS has been positively identihed in several cases (e.g. Frail et al. 1997). The 
subject of ISS is a large one (see review by Rickett 1990), and encompasses a range of 
observational phenomenology. Here we concern ourselves with estimating the importance 
of ISS-induced modulation of our flux density measurements as a function of time and 
frequency. For more detailed treatments of ISS see Goodman (1997) and Walker (1998). 

The character of the expected flux density fluctuations depend on whether the 
scattering takes place in the “strong” or “weak” regimes. We dehne a transition frequency 
z/q below which the density irregularities can induce strong phase variations and above it 
there is only weak scattering (WISS). Strong scattering can result in both slow, broad-band 
changes, due to refractive interstellar scintillation (RISS), and fast, narrow-band variations 
due to diffractive interstellar scintillation (DISS). The value of z/q along any given 
line-of-sight depends on the strength of the scattering and the total pathlength. The 
amount of scattering is parameterized by the scattering measure, SM, as modeled for the 
Galaxy by Taylor and Gordes (1993). For convenience we assume that the scattering along 
the line-of-sight takes place in a thin screen located at dgcr ~ {hz/2) sin|6|“^, where hz (— 1 
kpc) is the exponential scale height of the ionized gas layer (Reynolds 1989). 


Radio afterglow emission will be modulated by ISS provided that the source size 
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Og is close to one of the three characteristic angular scales (1) - (2) 

6d = and (3) 6r = 6fq{i^o/■, where 6p^ is the angle subtended by the first 

Fresnel zone at the scattering screen for v = z / q , i.e. 9fq = {c/2'KVQdscrY^‘^■ The modulation 
index mp (i.e. the r.m.s. fractional flux variation) is given by: 

for z/ > z/q and Og < 0^ (WISS), 

1 for z/ < z/q and Og < 9d (DISS), (1) 

(z/o/z^)“^'^'^^° for z/ < z/q and 9s < 9r (RISS). 

The expansion of the afterglow may cause 9g to eventually exceed one or more of the 
three characteristic angular scales that dehne the different scattering regimes. For ease of 
comparison in Appendix A we compute 9g for three different interstellar medium (ISM) 
models (a constant density ISM model, a wind model [WIND] and a jet model [JET]; See 
§3.2). Then the modulation will begin to “quench” as follows: 

{ {9w/9sY^^ for z/ > z/q and 9s > 9^ (WISS), 

{9d/9s) for u < uq and 9s > 9d (DISS), (2) 

{9r/9gY^^ for z/ < z/q and 9s > 9^ (RISS). 

For GRB991208 at (/,6) = (72°4,+42°6) we estimate SM = 2.26xl0“^kpc from 

the Taylor and Cordes (1993) model. Taking a typical scattering screen distance of dscr 
= 0.74 kpc we compute z/q = 4.8 GHz and 9fo = 4.2 pas (Walker 1998). The difficulties 
in obtaining accurate estimates of SM and dgcr likely limit the accuracy of the z/q and 9fq 
values to ± 30%. In Figure 2 we plot these values along with lines showing the different 
ISS regimes. On the right hand side of the hgure we compute a mo del-dependent timescale 
for the hreball to reach the angular size given on the left hand side. We have assumed that 
the emission from GRB 991208 originates from a jet-like outflow for tj^t < 2 days and used 
the relevant expression for 9s in Appendix A. At a redshift = 0.7055 (Dodonov et al. 
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1999, Djorgovski et al. 1999) the angular distance Da = Dl/{1 + ;2)^=0.45xl0^® cm, where 
Dl is the luminosity distance. The source size is Og = 5.2(ii^52/ni)^/®(trf/15)^/^ /ras, where 
E ^2 is the energy in the blast wave in units of 10^^ erg and rii is the density of the ambient 
medium in particles cm“^. Thus the observations in Tables 1 and 2 (and indicated as small 
crosses in Figure 2) he on either side of z/q and were taken during a time when dg expanded 
to eventually exceed 6p^. 

It is apparent in Figure 2 that most of the observations were taken in the weak 
scattering regime (z/ > z/q). Furthermore, for z/ > 15 GHz, dg > 9^ and therefore we expect 
that the observed variations will be dominated by source evolution and measurement error, 
not ISS. ISS starts to become comparable to the measurement error at 15 GHz where it is 
straightforward to show that mis = 0.08(trf/15)“’^/^^ (for Og > 9^). 

Accounting for ISS is particularly important at 8.46 GHz since the measurement 
error at this frequency is small (< 5%) while we expect m 8.46 = 0.45 for 9g < 9^ and for 
tci ^ 5 days m 8.46 = 0.24(G/15)“^/^^. Since the 8.46 GHz light curve is so well sampled, 
we can independently check on this prediction by subtracting out any linear trend and 
computing the scatter. For different assumed curves we hnd within the first two weeks of 
observations an r.m.s. scatter of 20-30%, suggesting that the strength of scattering along 
this line-of-sight may be slightly weaker than we have predicted. 

We also expect large variations in the flux density at z/ = 4.86 GHz but these are 
difficult to quantify, owning to the proximity of z/q (Walker 1998). For the purposes of 
model htting we adopt m 4,86 = 0.5 and for G ~ 10 days m 4.86 = 0.40(G/15)“’^'^^^. 

The line 9p = indicates the transition diameter that separates the 

region where the modulations are RISS dominated (above 9p) from that were they are DISS 
dominated (below 9t)- The observations at z/ = 1.43 GHz are clearly in the strong regime 
(z/ < z/q) where 9^ < 9g < 9r, and where RISS is expected to be signihcant with nip = 0.5 
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on refractive timescales of a day or more. Shorter timescale (foiss ~3 hrs) variations due to 
DISS are also expected but at a reduced level of 10-15%, suppressed by the hnite source size 
(i.e. 9s » Od) and bandwidth averaging. DISS, unlike WISS and RISS, is a narrow-band 
phenomenon with a decorrelation bandwidth Aud = < 100 MHz, the bandwidth 

of our observations. 


3.2. Synchrotron emission 

The radio/mm/optical to X-ray afterglow emission is believed to arise from the forward 
shock of a relativistic blast wave that propagates into the circumburst medium (see Piran 
2000, van Paradijs, Kouveliotou & Wijers 2000 for reviews). The electrons are assumed to 
be accelerated to a power-law distribution in energy, diV/dye cx for y > ym, and radiate 
synchrotron emission; where ye is the Lorentz factor of the electrons, ym is the minimum 
Lorentz factor, and p the power-law index of the electron energy distribution. Detailed 
calculations by Granot, Piran and Sari (2000a,b) and Granot and Sari (2002) have shown 
that to good approximation (to within a few percent) the synchrotron afterglow spectrum 
can be described by the following function 


X 

Where, Fm is the synchrotron peak flux density at the peak frequency z/m, and z/a is 
the synchrotron self-absorption frequency. The index n = —0.837 (Granot & Sari 2002) 
and controls the sharpness of the peak at z/^. We have assumed slow cooling, i.e. z/c, 

where z/^ is the synchrotron-cooling frequency. Gastro-Tirado et ah (2001) have measured 


F^(2)-/^[l-e-(^)-^^M^ 

\ 


nj'i 


5/3 


+ 


-{p-l)n/2 


Ijn 


(3) 
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a steeper optical slope f3 than Sagar et al. (2000) and have argued that this is evidence 
that the afterglow observations of GRB 991208 require the synchrotron cooling break to be 
located between radio and optical wavelengths. We model a possible synchrotron cooling 
break by a sharp spectral break with decrement A/? = 0.5 at z/ = Vc, and parameterization 

z/c = 

The evolution of the break frequencies z/a and z/m, z/c and the peak flux density 
follows from assumptions on the dynamics of the relativistic blast wave. It is common 
practice to assume three cases: i) propagation into a constant density ambient medium 
(ISM; Sari, Piran and Narayan 1998), ii) propagation into a wind-stratihed ambient medium 
(WIND; Chevalier and Li 1999, 2000), and iii) a blast-wave that is collimated in a jet (JET; 
Sari, Piran and Halpern 1999; note that the evolution in the case of a jet is rather insensitive 
to the density structure of the ambient medium; for details see Kumar and Panaitescu 
2000). The three cases have in common that the synchrotron parameters all evolve as power 
laws in time {i.e., F^a = , z/a = Cat”""*, z/m = and z/c = Cct”"'") but differ in 

the values of their decay indices and coefficients. In addition to assuming model-dependent 
decay indices, in this paper and in Paper I we solved for ap, a^., and Oc explicitly by 
model fitting the above synchrotron spectrum (FREE). The advantage of this approach is 
that it assumes only that the afterglow is due to synchrotron emission with a power-law 
distribution of electron energies, and it imposes no physical constraints on the geometry 
of the outflow, the microphysics of the shock or the density structure of the circumburst 
medium. 


3.3. Modeling the light curves 

We fit the function F(z/; Cp, Ca, Cm, Cc, ap, Oa, «m, «c,p) in Eqn. 3 to the observations 
presented in Paper I, this work (Tables 1 and 2) and the early time optical data of 
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Castro-Tirado et al. (2001). We minimize the of the fit by the Levenberg-Marquardt 
least-squares minimization scheme. We account for the effect of interstellar scintillation 
by adding in quadrature the predicted scatter due to ISS (see previous discussion) to the 
measurement error. We determine the 68.3 % (la) uncertainties in the parameters by 
plotting contours in chisquared where is raised by 2.3 as a function of the two relevant 
parameters and by leaving the remaining parameters free. The result of the best hts are 
shown in Figure 1 and presented in Table 3. 

In addition to this free-form model (FREE), we also hxed the decay-rate parameters 
(a^jCTa, «m, and Oc) at the predicted values for the constant density ISM, the WIND and 
the JET model and repeated the £t. The results of the JET fit are shown in Figure 3 
and presented in Table 3. As expected from the earlier attempts in Paper I, the standard 
WIND and constant density ISM fits were rejected as unsuitable (Xr > §)• These fits do not 
account well for the combination of the rapid decay at optical wavelengths and the more 
shallow decay at radio wavelengths. 


4. Discussion 

The principal challenge in modeling the afterglow of GRB 991208 is in simultaneously 
fitting for the decay of the radio (or = 1.07 ±0.09) and optical (oo = 2.2 ±0.1) light curves. 
The model that best represents the data is the free-form model (FREE) in which all the 
parameters in Eqn. 3 were solved for. The FREE model shown in Figure 1 and summarized 
in Table 3 differs from the original model presented in Paper I in several respects. Most 
notably, with the addition of the optical data from Castro-Tirado et al. (2001), it has been 
necessary to allow for a cooling break between the optical and millimeter bands (§3.2). 
This has the effect of flattening the value of p which improves the fit to the decay of the 
optical light curve. A cooling break in the spectrum also alters the temporal decline of the 
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synchrotron parameters but the values in Table 3 agree, within the errors, to those derived 
in Paper I. The FREE model resolves the differences between ctr and cTo by allowing for a 
relatively shallow decay of the peak flux density, i.e. Em oc and an unusually rapid 

evolution of the cooling frequency to lower frequencies, i.e. Vc oc This is generally 

much faster than any model yet proposed for GRB afterglows. For oc this implies 
F X R^=constant (where F is the bulk Lorentz factor of the flow, and B is the post-shock 
magnetic held). For any reasonable model both F and B are decreasing with time unless 
either the circumburst density or the fraction of shock energy in the B-£eld cb are increasing 
in time. In a future paper (S. Yost, in preparation) we plan to explore such modifications 
to the standard model for a number of well-studied afterglows. 

Both Li and Chevalier (2001) and Panaitescu & Kumar (2002) suggested a revision 
to the basic afterglow model in order to overcome the difficulties with jointly htting the 
radio and optical data for this burst. They proposed the existence of a two-component 
electron energy distribution with different slopes above and below some break frequency 
z/f,. Reasonable fits were obtained for either WIND and JET models provided that the 
radio (optical) emission was dominated by the low (high) energy electrons. Li & Chevalier 
(2001) also made predictions for the future evolution of the afterglow beyond the hrst two 
weeks, allowing further tests of their model. The measured ctr ~ 1.1 for the 8.46 CHz light 
curve is in good agreement with their predicted decay q;r = 1.25. The largest source of 
disagreement is in the low frequency behavior. In their WIND model the peak hux density 
at 1.43 CHz is predicted to be a factor of 10 times larger than what is measured. This is 
likely due to the relatively slow evolution of z/a where eta = 0.29, a value that is closer to the 
JET model (eta = 0.2) than the WIND model (eta = 0.6). 

In Paper I when the spectral fits were compared with the expectations for specific 
models, satisfactory agreement was found for the JET model. Nevertheless, there were a 
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number of discrepancies, notably the decline in the peak flux density Fm was too shallow, 
and the value for the slope of the electron energy distribution p = 2.52 was steeper than 
that implied by the decay of the optical light curves p ~ 2.2. It was proposed that these 
problems with the JET model could be reconciled if the jet was in transition to its fully 
asymptotic behavior. In our new JET model (see Fig. 3 and Table 3) the addition of a 
cooling break improves the overall quality of the £t considerably. Despite this agreement, 
there is a suggestion in Fig. 3 that the JET model is having some difficulties with the late 
time radio measurements. Most of the excess fhe JET model results from the fact 

that the late-time light curve (At >100 days) at 8.46 GHz is flat, q;r ~ 1.1, whereas the 
JET model predicts that the flux density will fall as F^, oc t~^, in agreement with cto ~ 2.2. 

There are two additional effects, not accounted for in our analysis, that could make 
«R < tto- In the first instance, GRB 991208 could have occurred in a radio-bright host 
galaxy {e.g., Berger, Kulkarni & Frail 2001). Sokolov et ah (2001) have used various optical 
methods to estimate the star formation rate (SFR) for the host galaxy of GRB 991208 and 
derive an extinction-corrected SFR> 100 Mq yr“^. Although there are large uncertainties 
in such estimates, at 2 : = 0.7 such elevated SFR’s should be detectable at centimeter 
wavelengths. Assuming a standard starburst spectral energy distribution (Yun & Garilli 
2002) we predict iG(1.43GHz) 110 pJy, significantly brighter than our measurements. 

We can therefore infer that the SFR in the host of GRB 991208 is probably no higher than a 
few tens of Mq yr“^. Deeper radio observations taken a year or more after the burst, when 
the emission from the afterglow is negligible, are needed to provide a more accurate estimate 
of the star formation rate. Another possible source of the flattening of the radio light curves 
could be the transition of the relativistic blast wave to non-relativistic expansion. Frail, 
Waxman, & Kulkarni (2000) modeled the decay of the radio light curves for GRB 970508 
at At ^ 100 days as a non-relativistic expansion. The expected decay rate in this phase of 
the evolution is q;nr = 3(p — l)/2 — 3/5 ~ 1.2 for p = 2.2. Thus, it is possible that the 



difference in decay rates is due to a transition to non-relativisitic expansion at t ~ 40 days. 
More sophisticated modeling is needed which explicitly includes the evolving dynamics of 
the blastwave {e.g., Harrison et al. 2001, Panaitescu & Kumar 2002) in order to test this 
hypothesis. 

In summary, with continued monitoring of the bright afterglow of GRB 991208, it has 
been possible to track the evolution of its radio light curve out to nearly 300 days after 
the burst. The power-law decline of the optical light curve at early times is considerably 
steeper than the radio light curve at late times. There are several different explanations 
for reconciling this disparate behavior, including modifying the standard relativistic model 
with non-standard electron distributions, and with circumburst densities or equipartition 
parameters that increase with time. However, the simplest explanation which is 
consistent with the data and requires no significant modifications, is that the blastwave of 
GRB 991208 entered a non-relativistic expansion phase several months after the burst. 
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Fig. 1.— The radio to optical lightcurves of GRB 991208. Data are from this work (Table 
1) and Paper I. Optical observations are from Castro-Tirado et al. (2001). Open symbols 
indicate the VLBA data (Table 2). The solid line indicates the best free-form fit to the light 
curves and the dotted lines are the predicted r.m.s. scatter due to interstellar scintillation. 
See §3 for more details. 
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Fig. 2.— Interstellar scattering toward GRB 991208 with the three regimes of scattering: 
weak (WISS), difractive (DISS) and refractive (RISS). The solid lines indicate the frequency 
dependence of each of three characteristic angular scales, joined at the transition frequency 
uq and the Fresnel scale Opo- On the right vertical axis we compute a model-dependent 
timescale, corresponding to the angular radius reached by a collimated hreball given on the 
left vertical axis. The small crosses indicate the time and frequency of the measurements in 
Tables 1 and 2. See 53 for more details. 
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Fig. 3.— The radio to optical lightcurves of GRB 991208. The solid line indicates the best 
£t to the light curves for the JET model and the dotted lines are the predicted r.m.s. scatter 
due to interstellar scintillation. See 53 for more details. 
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Date 

(UT) 

tstart - tend 

(hh:mm:ss) 

tint 

(hh:mm:ss) 

At 

(days) 

F8.46 

(/4Jy) 

F 4.86 

(liJy) 

F 1.43 

(liJy) 

1999 Dec 23.80 

19:03:30 - 19:17:30 

11:20 

15.61 

782 ± 40 



1999 Dec 26.78 

18:25:40 - 18:47:20 

16:20 

18.59 

693 ± 38 



1999 Dec 28.98 

23:14:50 - 23:38:10 

35:30 

20.79 

464 ± 44 



1999 Dec 29.49 

11:17:10 - 12:11:50 

22:00 

21.30 

679 ± 30 



1999 Dec 29.49 

11:17:10 - 12:11:50 

18:40 

21.30 


501 ± 42 


1999 Dec 29.49 

11:29:20 - 12:41:17 

36:20 

21.30 



79 ± 41 

2000 Jan 5.94 

22:25:00 - 22:37:10 

9:30 

28.75 

646 ± 58 



2000 Jan 9.89 

21:04:50 - 21:48:00 

19:00 

32.70 

337 ± 34 



2000 Jan 9.89 

20:58:00 - 21:37:50 

14:25 

32.70 



61 ± 53 

2000 Jan 13.79 

18:13:10 - 19:33:00 

31:30 

36.60 

574 ± 25 



2000 Jan 13.79 

18:04:55 - 19:41:05 

30:35 

36.60 



56 ± 40 

2000 Jan 17.85 

19:33:40 - 21:17:30 

43:30 

40.66 

400 ± 23 



2000 Jan 17.85 

19:29:40 - 21:00:10 

26:30 

40.66 



73 ± 42 

2000 Jan 21.92 

21:17:30 - 22:36:20 

31:40 

44.73 

458 ± 31 



2000 Jan 21.92 

21:11:50 - 22:43:50 

28:30 

44.73 



98 ± 59 

2000 Jan 24.81 

19:03:25 - 19:52:00 

36:50 

47.62 


212 ± 29 


2000 Jan 25.85 

18:32:10 - 22:25:20 

1:00:10 

48.66 

442 ± 23 



2000 Jan 25.85 

18:28:00 - 22:19:40 

1:02:55 

48.66 



54 ± 37 

2000 Jan 30.85 

18:42:30 - 22:22:30 

1:20:50 

53.66 

290 ± 21 



2000 Jan 30.85 

18:38:10 - 22:13:20 

1:02:55 

53.66 



43 ± 41 

2000 Feb 4.80 

18:22:50 - 19:59:20 

38:00 

58.61 

239 ± 24 



2000 Feb 4.80 

18:18:30 - 20:06:30 

31:45 

58.61 



82 ± 43 

2000 Feb 10.78 

17:59:10 - 19:18:00 

31:30 

64.59 

232 ± 27 



2000 Feb 10.78 

17:55:00 - 19:26:40 

26:55 

64.59 



60 ± 49 

2000 Mar 13.63 

14:42:50 - 15:37:20 

40:50 

96.44 

172 ± 33 



2000 Mar 27.57 

13:34:20 - 13:48:50 

11:50 

110.38 

177 ± 47 



2000 Apr 13.67 

15:43:50 - 16:36:20 

38:50 

127.48 

123 ± 26 



2000 May 4.62 

14:27:40 - 16:20:50 

2:01:10 

148.43 

106 ± 23 



2000 May 7.64 

14:30:50 - 16:24:50 

1:24:20 

151.45 

152 ± 25 



2000 Jun 6.45 

10:31:10 - 11:09:30 

24:40 

181.26 

38 ± 40 



2000 Jun 14.40 

09:24:10 - 10:01:10 

27:00 

189.21 

29 ± 30 



2000 Aug 27.76 

22:15:30 - 23:08:10 

42:00 

263.57 

31 ± 27 



2000 Sep 10.85 

20:28:10 - 21:19:10 

01:25:20 

277.66 

44 ± 21 



2000 Sep 16.78 

18:37:00 - 21:25:30 

02:28:50 

283.59 

20 ± 19 



2000 Sep 24.77 

18:35:30 - 20:53:40 

02:02:03 

291.58 

51 ± 15 



2000 Sep 27.10 

02:19:10 - 03:40:30 

01:05:30 

293.91 

48 ± 20 




tstart and end 
tint, the time 
Observations 

taken before Dec 22 1999 can be fonnd in Paper I. 


Table 1: Snmmary of VLA observations of GRB 991208. The date, start time 
time fend of the observations are given in UT. Also, the total time on sonrce 
passed since the event At and the flnxes at 8.46, 4.86 and 1.43 GHz are given. 
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Date 

(UT) 

At 

(days) 

Fs.ss 

(z^Jy) 

F4.85 

(/iJy) 

1999 Dec 14.73 

6.54 

967 ± 55 

448 ± 67 

1999 Dec 18.72 

10.53 

1008 ± 60 

888 ± 74 

1999 Dec 26.70 

18.51 

430 ± 75 

431 ± 92 

1999 Jan 4.69 

27.50 

474 ± 56 

440 ± 64 


Table 2: Summary of VLBA observations of GRB 991208: the date, the time passed since 
the event At and the fluxes at 8.35, and 4.85 GHz. 


Table 3. Best Fit Model Parameters for GRB 991208 : the value of p, the peak flux 
density Fm = Cp x , the synchrotron self-absorption frequency z/a = Ga x the 
peak frequency z/m = Gm x the cooling frequency z/c = Gc x the reduced 
chisquared Xr the degrees of freedom (d.o.f.) of the £t. 


p Cf Cl, Cm Cc ap Qfa am ttc Xr d-O-f Hote 

(mJy) (GHz) (THz) (THz) 


2 14+0-16 

7 7+15.7 
‘ •' -2.3 


0 34+0-29 

^•o^_0.30 

2670^1^ 

0 

^"^^-0.74 

0 

'^•^^-0.17 

1 q+i-8 

i.y_l 2 

2.081+? 

1.52 

104 

FREE 

2 ^f^+O-io 

17.8+J 

6.3+;° 

3Q 4+00.1 
oy.^_3 10 

3.0+;? 

1.0 

0.2 

2.0 

0.0 

1.89 

108 

JET 


Note. — Fixed parameters can be recognized by the fact that they have no quoted uncertainties. 
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A. Observed Angular Size 

The relation between observed time, T, real time t, the radius R and the angle 
/i = cos 9 is given by 


T = t — nR/c. (Al) 

R and t are related by some hydrodynamic solution given by the dynamics of the shock, 
so that in a given observed time T, i? is a function of jj, only. Shock dynamics are usually 
described by a powerlaw, 7 R Therefore, 

1 

^ 4(m + 1 ) 7 ^ 

Substituting into equation Al gives 



CT / R = 1— yU + 


4(m + 1 ) 7 ^ 

Following Sari (1998) we dehne 7 ^ to be the Lorentz factor of the fluid behind the shock on 
the line of sight and Rl to be the shock radius on the line of sight, both at time T. We get 

Rl 


4(m + 1 ) 7 ! 

Substituting this in the previous expression, we have 


(A 2 ) 


1 = ;^4(m + 1)7! 


1 — P- + 


4(m+ 1 ) 71 , \Rl 


R 


Which can be rewritten as 

1 — /i = 


4(m + 1 ) 7 ! 


R 


-1 


R 

Tl 


The perpendicular size, in which we are interested, is given by 

R_l = RsinO = R^/l — = \/2R^l — fi. 
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It is therefore maximal when R^{1 — fi) is maximal. 


^ V ( 4 -) 


m+2 


and the maximum is obtained for 


R 


Rl \Rl 


l/(m+l) 


Rl \m + 2^ 

The maximal perpendicular size from which radiation arrives is therefore given by 


R 


Ry = M-ET 

7 l \RlJ \ 4(m+l) [\Rl 


-1 


R 

Tl 


RlHl 


l/2(m+l) 


^/2{rnT^ \m + 2^ 

We now only need to get the coefficient for the relation between Rl and T to get a closed 
expression. Roughly this is given by i? ~ R?pc^'y‘^ where p is the density just in front of the 
shock. The more exact result is given by Blandford and Mckee (1976; their Eq. 69) as 

167r 


E = 


5 + 4m 


R^pc^j'^. 


(A3) 


If we write p = ARJ^ this, together with equation A2, translates to 




Rl 


5 + 4m 4(m + l)cT 


Rl/il = 2 a/ (m + l)cT 


Rl = 


(5 + 4m)(m + 


47rcA 


Therefore: 


R±,max 


2(m + l)cT 
m + 2 


(5 + 4m)(m + l)Ti7 
47rcA 


(5 + 4m) (m + l)Ti7 
47 rcA(m + 2) 


l/2(m+l) 


l/2(m+l) 


For a wind we have m = 1 so 


Rl 


max,WIND 


4cT 


ATE 

2t[cA 


1/4 


= 1.1 X lO^^cm 


1 + ;^ 


-3/4 


Rj/A:R/T/15days)3// 
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where E = 10^^ii^52, we have corrected for redshift and have used A = 5 x g cm“ 

as in Chevalier and Li (2000). Or in angular size 

0WIND = 2.2/ias D^28EU"A;^/\T/15 de.ysf\ 

where L*a ,28 is the angular distance in units of 10^® cm. For ISM models we have m = 3, 
and 

/I _L \ “5/® 

OisM = 2.8/ias {-^—J days)^/®, 

where ni is the density of the ISM in particles cm“®. We can now obtain an estimate 
of the angular diameter for a jet as follows. During sideways expansion of the jet, 
the radius R of the hreball is practically constant, 7 oc (Sari, Piran & Halpern 

1999) and so i?_L,max,jET cc At the time of the jet break tjEx we thus roughly have 

max,JET ~ Rt, max,ISM (Oet) X (f/tjEx)^^^, i-e. 

/ 1 -L. \ “5/® 

6'jet = 1.7/jas {-^—j -0)4,28-^52*% ^'^*(OET/8hr)^/®(T/15days)^/^. 



